Equivalence of driven and ageing fluctuation-dissipation relations in the trap model 
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We study the non-equilibrium version of the fluctuation dissipation (FD) relation in the glass phase 
of a trap model that is driven into a non-equilibrium steady state by external shear. This extends 
our recent study of ageing FD relations in the same model, where we found limiting, observable 
independent FD relations for "neutral" observables that are uncorrelated with the system's average 
energy. In this work, for such neutral observables, we find the FD relation for a stationary weakly 
driven system to be the same, to within small corrections, as for an infinitely aged system. We 
analyse the robustness of this correspondence with respect to non-neutrality of the observable, and 
with respect to changes in the driving mechanism. 
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I. INTRODUCTION 

Glasses relax very slowly at low temperatures. They 
therefore remain out of equilibrium for long times and 
exhibit ageing 0: the time-scale for response to an ex- 
ternal perturbation (or for the decay of correlations) in- 
creases with the "waiting time" since the system was 
quenched to the low temperature, and thus eventually 
far exceeds the experimental time-scale. Time transla- 
tional invariance (TTI) is lost. As a result of this dy- 
namical sluggishness, glassy systems are highly suscepti- 
ble to external driving, even when the driving rate 7 is 
small. One example of 7 is shear rate in a rheological 
system. Typically, steady driving interrupts ageing and 
restores a non-equilibrium steady (TTI) state in which 
the time-scale defined by the inverse driving rate plays 
a role analogous to the waiting time of the ageing 
regime |l|,||, p. 

Ageing and driven glasses in general violate the equi- 
librium fluctuation dissipation theorem (FDT) |2l]. Con- 
sider the autocorrelation function for a generic ob- 
servable m, defined as C{t,t^) — {m{t)m{t^)) — 
{m{t)) {m{t^)). The associated step response function 
x(i,iw) = dt' R{t,t') tells us how m responds to a 

in its conjugate field h. In 
- tw) by TTI (similarly for 
:Xit - iw) = R{t - iw) = 

<5(m(t)> 



small step h{t) — hQ{t — t„ 
equilibrium, C{t,t^) ~ C{t 
x), and the FDT reads —-. 



T dU 



■C{t - iw) where R{t — iw, ^w) 



is the 



h=0 



impulse response function and T is the thermodynamic 
temperature. (We set /cb = 1-) A parametric 'FD plot' 
of X vs. C is thus a straight line of slope — l/T. 

Out of equilibrium, violation of FDT is measured by 
an FD ratio, X{t, t^), defined through |l| 



_d_ 



X(t,t^) d 



C{t, u 



(1) 



long times, indicating strongly non-equilibrium behavior 
even though one-time observables such as entropy may 
have settled to essentially stationary values. Similarly, 
driven glassy systems can violate FDT even in the limit 
of weak driving. 

Remarkably, the FD ratio in several ageing mean field 
models |l2| assumes a special form at long times. 
Taking tyj ^ 00 and i — s- 00 at constant C = C{t,t^), 
X{t, tyj) ^ becomes a (non-trivial) function of the 

single argument C . If the equal-time correlator C(t, t) 
also approaches a constant Cq for t ^ cxd, it follows that 



1 /■'^o 

X(i,tw) = 7^ / dCX 
^ Jc(t,t^) 



:{C). 



(2) 



Graphically, this limiting non-equilibrium FD relation 
is obtained by plotting x ^s. C for increasingly large 
times. From the slope —X{C)/T of the limit plot an 



In ageing systems, violation {X ^ 1) can persist even at 



effective temperature |14| Tcs{C) = T/X{C) can be de- 
fined. Throughout this paper, we absorb the factor T into 
the response function so that an equilibrium FD plot has 
slope —1. 

In the most general ageing scenario, a system displays 
dynamics on several characteristic time scales, each with 
its own functional dependence on t^j. If these time scales 
become infinitely separated as tw ^ 00, they form a set of 
distinct 'time sectors'. In mean field, Tcf[{C) is constant 
within each sector p^ ] , and independent of the observable 
TO used to construct the FD plot. These properties have 
also been observed in some lower dimensional (non mean 
field) systems 

Gugliandolo et al. lIT^ proposed that an equivalent lim- 
iting non-equilibrium FD relation should hold in slowly 
driven glassy systems (7 — + 0) and that the corresponding 
effective temperatures TofflC', 7 — ^ 0) and Tofi(C, tw — > 
00) should coincide. Although this is believed to apply 
widely among driven glasses, the evidence supporting it, 
to date, is limited to the two detailed studies of Berthier, 
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Barrat et al., in mean field g and in simulations of 
sheared Lennard Jones particles, initially in Ref. ||] 
and laterj^with a study of observable independence, in 
Refs H, §. We note that FDT has also been studied 
in a driven phase separating model (T < Tc where Tc 
is the critical temperature) with N non-conserved order 
parameters in the large N limit although here the 
drive does not interrupt ageing. In this study, in the 
timesector with stationary dynamics, the FD plot is of 
trivial equilibrium form for 7 = but a non-trivial curve 
for 7 > 0; in the ageing sector, the response function is 
constant, giving a flat FD plot. 

In this paper, therefore, we study FDT in the driven 
regime of Bouchaud's trap model ^ |2^, for 
which correlation and response functions can be calcu- 
lated exact ly. This will allow a comparison with our re- 
cent study pSl of the ageing trap model, where we found 
limiting FD relations that are observable independent for 
"neutral observables" that are uncorrelated with the sys- 
tem's energy. This is consistent with the mean field work, 
for which observables are usually defined in terms of ran- 
dom couplings, uncorrelated with the average energy. (In 
coarsening models, similar arguments have been used to 
exclude observables correlated with the order parame- 
ter jl].) Surprisingly, however, we found the FD plot to 
be a continuous curve even though the model has just 
one time sector, with relaxation times 0(iw)- Although 
this finding is apparently at odds with the mean field 
predictions, it is likely to result from the fact that the 
trap model has a broad distribution of relaxation times 
(all within its single time sector). We return to this point 
in the conclusion. 

In what follows, our central result will be that the same 
non-trivial FD relation is found (to within logarithmic 
corrections), even when the trap model is weakly driven 
according to the mechanism proposed in Ref. Al- 
though the curvature of the FD plot obviously excludes 
a constant effective temperature, our finding is still con- 
sistent with Cugliandolo et aZ.'s predictions, in so far as 
the relationship between correlation and response is the 
same for aged and weakly driven glasses. This finding is 
non-trivial since the shapes of the relaxation spectra for 
the ageing and weakly driven trap models differ strongly 
from each other. 

We start (Sec. ||) by defining the trap model and sum- 
marising the ageing FD predictions of Ref. . We then 
(Sec. |l|) derive exact expressions for the autocorrela- 
tion and response functions for an arbitrary observable 
m in the steadily driven regime. Using these, we calculate 
the limiting driven FD relation, for neutral observables 
that are uncorrelated with the average energy. We show 
that this relation is the same (to within logarithmic cor- 
rections) as its ageing counterpart (Sec. [V). We then 
(Sec. ^ consider robustness of this correspondence with 
respect to (i) changes in the driving mechanism and (ii) 
non-neutrality of the observable, before concluding. 



II. THE TRAP MODEL; DRIVING 
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z=E/T-Iog(t^) or z=E/T+log(y) 

FIG. 1: Scaled energy distributions for (i) the ageing model 
at = 10^, 10", 10^, 10'^ (dashed lines) and (ii) the driven 
model for 7 = 10'^, 10"", 10-^ 10"® (solid lines; different 7 
values indistinguishable). The temperature T = 0.3. 

The trap model comprises an ensemble of uncoupled 
particles exploring a spatially unstructured landscape of 
(free) energy traps by thermal activation. The tops of 
the traps are at a common energy level and their depths 
E have a 'prior' distribution p(E) {E > 0). A particle 
in a trap of depth E escapes on a time scale t{E) = 
To exp{E/T) and hops into another trap, the depth of 
which is drawn at random from p{E). The probability, 
P{E,t), of finding a randomly chosen particle in a trap 
of depth E at time t thus obeys 

dtPiE, t) = ^T-\E)P{E, t) + Y{t)p{E) (3) 

in which the first (second) term on the RHS repre- 
sents hops out of (into) traps of depth E, and Y{t) = 



{r-\E)) 



is the average hopping rate. 



For the specific choice of prior distribution p{E) ~ 
exp(— i?/Tg), the model shows a glass transition at a 
temperature Tg. This is seen as follows. At a temper- 
ature T, the equilibrium state (if it exists) is Peq{E) oc 
t{E)p{E) oc exp (E'/T) exp (— i5/Tg). For temperatures 
T < Tg this is unnormalizable, and cannot exist; the 
lifetime averaged over the prior, (t)^, is infinite. Fol- 
lowing a quench to T < Tg, the system never reaches a 
steady state, but instead ages. This can be seen from 
the time evolution of P{E,tv,), which can be obtained 
exactly from Eqn. ^ (with Y{tyj) determined self consis- 
tently by enforcing normalisation of P{E, t^) |2^). At 
large times 00 a scaling limit is reached in which 

Pij,t^) = \T /T{E)\P[E,t^) is concentrated entirely on 
traps of lifetime r = 0{t^): the scaling distribution 
P{z) = TP{E) where z = E/T - log(iw) is shown in 
Fig. The model thus has just one characteristic time 
scale, which grows linearly with the age iw (In contrast, 
for T > Ta all relaxation processes occur on time scales 
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0(to).) In what follows, we rescale all energies such that 
Tg = 1, and times so that Tq = 1. 

Driving was first incorporated into the model in or- 
der to study the rheology of "soft glassy materials" . Al- 
though we are not directly interested in rheology here, we 
use the same driving rules which are defined as follows. 
Each particle is assigned its own local elastic "strain" 
with a corresponding "stress" kl. Each time the particle 
hops, I is set to zero. Between hops, / = 7 where 7 is 
the rate of external driving ("straining"). A particle in 
a trap of depth E strained by I sees a reduced effective 
energy barrier E 



^kP, so that 



DtPiE,l,t) 



-\E)e''^'/^^ P + Y{t)p{E)S{l). (4) 



In this equation, T^^{E)e^^ ^"^"^ is the strain-enhanced 
counterpart of the "bare" activation rate T^^iE) defined 
above for the undriven model, and Dt is the convected 
derivative Dt — dt+'^di. Eqn. ^ (integrated on I) reduces 
to Eqn. ^ for 7 = 0, as required. In the following we 
rescale I such that k — \. 

For steady driving (7 = const.) ageing is interrupted, 
and a TTI steady state is restored. To write down the 
steady state distribution, define S{li,l2,E) as the prob- 
ability for a particle that starts off with strain li and in 
a trap of depth E not to hop until its strain has reached 
l2- From Eqn. this is 



S{lij2, E) — exp 



-{E)i 



dl exp 



2T 



(5) 



In terms of this quantity, the steady state distribution of 
Eqn. I can then be written as 



Poo{E,l) 



■piE)S{0,l,E). 



(6) 



Here 



7 



is the steady state average hopping rate. 



which can be determined self consistently by enforcing 
normalisation of P^{E,l). (This is most conveniently 
done by changing variables to t = exp(i?/T) which gives 



P{r) 



.-l-T 



) In the limit 7-^0, the energy distribu- 



tion Poo{E) — J dl Poo{E, I) approaches a scaling limit in 
which all relaxation times are 0(1/7). The scaling distri- 
bution Poo{z) = TPoo{E) where z = E/T + log{'^) differs 
strongly from its ageing counterpart P{z), as shown in 
Fig. I 



III. CORRELATION AND RESPONSE 

FDT can be studied by assigning to each trap, in addi- 
tion to its depth E, a value for an (arbitrary) observable 
m The trap population is then characterized by 

the joint prior distribution a{m\E)p{E), where a{'m\E) 
is the distribution of m across traps of given fixed energy 
E. The dynamics then obey 



DtP{E,m,l,t) 



-\E, m)e'=''/2T p ^ Y{t) p{E)6{l)a{m\E) 



(7) 



where the activation times are modified by a small field 
h conjugate to m as T{E^m) = t{E) exp (mh/T). This 
particular form of t{E, m) is one of several possible 
choices that all maintain detailed balance under zero- 
driving conditions |2^. We adopt it because, in the 
spirit of the unperturbed model {h = 0), it ensures that 
the jump rate between any two states depends only on 
the initial state, and not the final one. 

In Ref. we derived exact expressions for the 

two-time autocorrelation and step response functions, 
C{t,t„) and x(tjtw) in the ageing regime (t„ 00, t 
00 at fixed t/t^) of the undriven model (7 = 0), following 
a quench into the glass phase at time — from an ini- 
tially infinite temperature. Each comprises two compo- 
nents that depend separately upon the functional forms 
of the mean, rn{E), and variance, A'^{E), of the distri- 
bution a{m\E). For the purposes of this paper we are in- 
terested only in observables with rn{E) = and (within 
these) mainly the neutral observable for which the vari- 
ance is uncorrelated with energy, A^(i?) = const. = 1. In 



this case, in the simultaneous limit ty, ^ 00 with t — > cx), 
C{t,tvj) and x(^?iw) depend on time only through the 
scaling variable {t — t^)/t^ as shown in Fig. ||a,b. The 
corresponding FD plot is shown in Fig. |^. 

In the steadily driven regime, TTI is restored: C and 
X do not depend explicitly upon the waiting time but 
only on the measurement interval t — iw, so we set iw = 
without loss of generality. For observables with rn{E) — 
0, the autocorrelation function is exactly 



dl 



7* 



dE A^{E)Poo{E,l). 



(8) 



This can be understood as follows. When any parti- 
cle hops, its new value of m is uncorrelated with the 
old one. At time t, therefore, only those particles that 
have not hopped since t = Q can contribute to the cor- 
relator, with weight A^(£'). The fraction of such par- 
ticles which had strain I and trap depth E at time 
t = is P^{E, l)S{l, I + 7t, E). From Eqn. | this equals 
Poo{E,l + 'jt) and integration on / and E gives the re- 
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suit, Eqn. ^. [Alternatively, Eqn. |8| can be understood 
by recalling the driving dynamics: upon any hop, each 
particle resets its local strain to zero; between hops, the 
local strain afhnely follows the applied one. Therefore, 
the particles that have not hopped since t = are just 



those that have strains / > 7*.] 

The corresponding switch-on response function has 
contributions from both hopped and unhopped particles. 
The contribution from hopped particles can be expressed 
as an integral over the last hop-time of each particle, t': 



Xhop{t,j) ^ dh\h=o I dm dE I dt' Y{h,t')ma{m\E)p{E) exp 
) Jo Jo 



7t(£', m) 



7(t-*') 



dsexp — - 



(9) 



r 



In this expression, Y{h,t) — Yoo + 0{h) is the fraction 
of particles that last hopped at time t' . Of these, a pro- 
portion a{m\E)p{E) chose energy E and magnetization 
m. The subsequent survival probability over the inter- 
val t' = ... t is encoded by the exponential factor. In 
principle, the differentiation on h has two contributions: 



J 



one from the factor T{E,m) in the exponential, and an- 
other from Y{h,t). However the second gives zero, since 
ma{m\E) — for the zero- mean variables consid- 
ered here. Adding the contribution from particles that 
have not hopped since i = 0, and doing some manipula- 
tion, we find finally the exact result 



dE 



r 



Poo{E,l + jt)] 1^ ds exp 



(10) 



(into which we have absorbed a factor T, as described 
above). 



IV. FD PLOTS 

Using the exact expressions of Eqns. || and |l^, we cal- 
culated C(i, 7) and x(t, 7) numerically for the neutral 
observable, A^(i?) = 1. As an independent check, we 
calculated each quantity by direct simulation, using a 
waiting time Monte Carlo technique. The results are 
shown in Fig. |a,b. In the limit 7 ^ 0, i ^ at fixed 
jt, C(t,'j) and x(t,7) depend on 7 and t only through 
the scaling variable 7t. Analytically, this can be seen for 
the correlator by substituting Eqn. ^ into Eqn. |^ and in- 
tegrating on dE by changing variables to r, as described 
above. The 7 dependence from the integral exactly can- 
cels the prefactor ^00/7 so that the only dependence on 
7 and t appears through the scaling variable jt in the 
limit of the integral on L A similar argument applies to 
the response function. 

The scaling functions C{'jt) and xiit) both differ 
strongly from their ageing counterparts (compare the 
solid and dashed lines in Fig. ||a,b). This is to be ex- 
pected, due to the obvious difference between the (scaled) 
energy distributions of the driven and undriven models. 
Remarkably, however, the driven FD relation x(C') is 
strikingly similar to its undriven counterpart (Fig. 0c). 
Both start with a slope ~X{C = 1) = x'(C =!) = -! 



(thus reproducing the equilibrium FD form in this limit) 
and finish with a slope x'(C' = 0) = at an intercept 
x{C = 0) = T. (We have confirmed these features ana- 
lytically as well as numerically.) Even between these lim- 
its, these is little discernible difference between the ageing 
and driven FD plots. This non-trivial result is consistent 
with the predictions of Cugliandolo et al, that the rela- 
tionship between correlation and response should be the 
same in weakly driven and old ageing glassy systems. 



In the inset of Fig. we show an expanded region 
of the main FD plot. Despite the striking similarity 
of the ageing and driven FD relations, our numerics do 
nonetheless suggest a small discrepancy. To investigate 
this further, we examined the behaviour of X = — x'(C) 
in the limit C ^ 0. Setting (t — t^)/t^ = u for 
the ageing model, we found C ^ u'"^ and X ^ u^^, 
hence X ~ C^/^ as u ^ 00. Setting -^t = v for 



the driven model, we found C 



exp{—v /2) and 



X ~ exp[-z;V2T], hence C ~ X^[\og{l/ X)]^^-^^^ as 
V —^ 00. Therefore, the driven and ageing FD plots are 
indeed equivalent in the limit C 0, but only to within 
minor logarithmic corrections. This explains the slight 
discrepancy seen in our numerical data. 
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FIG. 2: a) Correlator and b) response vs. scaled time for the 
neutral observable A'^{E) = 1 in the driven model (solid lines) 
and ageing model (dashed lines), calculated from the exact 
analytical expressions in the text and in Ref. [|l5|l. For the 



ageing case, waiting times = 10"^, 10*, 10^, 10^ are shown 
(but are indistinguishable from each other); for the driven 
case shear rates 7 = 10~^, 10~*, 10~^, 10~® are shown (also 
indistinguishable). As an independent check, we also show 
the driven correlator and response at 7 = 10~^ calculated 
from waiting time Monte Carlo simulation. For C this is 
indistinguishable from the exact results for C in a); for x it 
appears as the jagged line in b). (c) FD plots of correlator 
vs. response for the driven case (solid lines) and ageing case 
(dashed lines) constructed from the exact results of (a,b); the 
inset is zoomed on a small region of the main plot. For the 
driven case (solid lines), driving rate decreases downwards at 
fixed C. The temperature T = 0.3. 



V. ROBUSTNESS UNDER CHANGE OF 
DRIVING MECHANISM; NON-NEUTRAL 
OBSERVABLES 

We now investigate the robustness of the equivalence 
between driven and ageing FD relations with respect to 
(i) non-neutrahty of the observable m, and (ii) changes 
in driving mechanism. We start with (i), considering ob- 
servables for which A^(£') — exp{nE/T) (which defines 
n), though still with rn{E) = 0. In this case, the overall 
amplitude (initial value) of the correlator depends explic- 
itly on the waiting time (or driving rate) , even in the age- 
ing (or weakly driven) limit. In order to obtain a limiting 
FD plot, the correlator and response must be normalised 
by the initial value of the correlator. In the driven regime, 
therefore, we now plot x(i,7) = 7)/C'(0i 7) versus 
C{t,j) — C(i, 7)/C(0, 7), with t as the plotting param- 
eter. In the ageing case, for these non-neutral observ- 
ables, care must be taken in constructing the FD plot. 
For neutral observables, the usual prescription is to plot 
x(i,iw) versus C{t,t^) with t as the plotting parame- 
ter. For non-neutral observables, however, we see from 
Eqn. |l| that the slope of the FD plot is only guaranteed 
to coincide with X if we plot x{t,t^) = x(t, iw)/x(^j 
versus C{t,t^) = C{t,t^)/C{t,t) with as the plotting 
parameter. (This coincides with the usual prescription 
for neutral observables, as required.) These normalised 



FD plots are shown for n = 0.2 in Fig. ^, and are seen 
to differ strongly from each other: equivalence of ageing 
and driven FD relations does not hold in the trap model 
for non-neutral observables. 



O.lr 




FIG. 3: FD plots for (i) the ageing model with t = 10^, 10^ 
(dashed lines; different t values indistinguishable) and (ii) the 
driven model with 7 = 10"^, 10"", 10"^, 10"*^ (solid lines; 
different 7 values indistinguishable) for a non-neutral observ- 
able with m{E) = 0, A^{E) = exp(n£'/T) for n = 0.2. The 
temperature T = 0.3. 




FIG. 4: FD plots for the neutral observable with m(E) = 0, 
A^(£;) = 1 for (i) the ageing model with t = 10^ 10*, 10^ lO'' 
(dashed lines; different t values indistinguishable) and (ii) a 
driven model with the alternative driving dynamics of Eqn. ^ 
for 7 = 10-^ 10"*, 10-^ 10"'* (dot-dashed lines). The tem- 
perature T = 0.3. 



Finally, we consider robustness of the equivalence of 
ageing and driven FDT with respect to a change in the 
driving mechanism. To do this, we consider a differ- 
ent driving mechanism which adds an additional hopping 
process to the trap model whose rate 7 is independent of 
trap depth. This just has the effect of normalizing each 
hopping rate according to l/r — > I/T-I-7. The dynamics 
are now 



dtP{E,m,t) = - 



1 



t(E, m) 



+ 7 



P + Y{t) p{E)a{m\E). 

(11) 



6 



with a steady state energy distribution given by The calculation of correlation and response functions is 

now trivial since all survival probabilities are simple ex- 
p fTp\ _ ^ooP{^) /-,r,N ponentials; one finds 



POO 

Jo 



and 



X(t,7)=/ dEA'{E)P^{E) V (l-exp 
Jo l + -fT{E) [ 



r{E) 



+ i]t 



(13) 



(14) 



The FD plot for the neutral observable A^(i?) = 1 is 
given in Fig. ^, and seen to differ strongly from the cor- 
responding FD plot in the undriven model. Hence, the 
equivalence of the ageing and driven FD relations is not 
preserved for this change of driving mechanism. How- 
ever, this is likely to be a consequence of the fact that 
this second choice of driving mechanism does not in fact 
violate detailed balance, but merely renormalises all the 
jump rates. 

VI. SUMMARY AND CONCLUSION 

In this paper we studied the non-equilibrium FDT in 
the glass phase of Bouchaud's trap model []6|, pO| p2|, p3|, 
extended to incorporate the non-linear driving mecha- 
nism of Ref. |2^. After deriving exact expressions for 
the correlation and response functions of a generic ob- 
servable, m, we compared the FD relation for a system 
driven steadily at rate 7 — > with that for an ageing sys- 
tem at waiting time 1/7. For "neutral" observables that 
are uncorrelated with the system's average energy, the 
driven and ageing FD relations are the same, to within 
minor logarithmic corrections. This correspondence does 
not apply to non-neutral observables. Finally, we consid- 
ered an alternative driving mechanism that renormalises 
all the hopping rates according to 1/t ^ I/T + 7. In this 
case, the ageing and driven FD relations differ strongly, 
even for neutral observables. Although this is appar- 
ently at odds with our central result, it is likely to re- 
sult from the fact that this trivial driving mechanism 
does not violate detailed balance. Further research is 



certainly needed, however, to understand whether other 
conditions are required on driving mechanisms in order 
to get steady-state behaviour related to that of undriven 
ageing systems. 

We return finally to address the fact that the FD rela- 
tions are rounded in this model, thus excluding a sin- 
gle effective temperature within the ageing or driven 
time sectors (with relaxation times 0(tw) or 0(7) re- 
spectively). Similarly 'rounded' FDT plots have recently 
been found in coarsening models at criticality | [l6t ; the 
limiting value —Xq^ of the slope for C ^ was there 
shown to be a universal amplitude ratio. It is possible 
that at least this X^o could define a sensible Tea, and 
in fact both our limiting FDT plots (for the first driving 
mechanism) share a common value X^o — 0. 

In conclusion, the FD relation between correlation and 
response is the same, to within logarithmic corrections, 
in the ageing and driven trap models. We suggest that 
the limiting value — Aoc of the slope for C — > could be 
used to define an effective temperature. Further work is 
needed to delineate more fully the class of finite dimen- 
sional driven glassy models that exhibit this behaviour. 

Note added: we have recently become aware that Lu- 
dovic Berthier has studied FDT in the driven EA model, 
and found FD plots the same as in the ageing model, 
with agreement between two different observables (un- 
published) . 
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